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Abstract 

On ■ 

We prove a convergence theorem for a sequence of super-Brownian motions moving among hard 
Poissonian obstacles, when the intensity of the obstacles grows to infinity but their diameters shrink 
to zero in an appropriate manner. The superprocesses are shown to converge in probability for the law 
P of the obstacles, and P-almost surely for a subsequence, towards a superprocess with underlying 
spatial motion given by Brownian motion and (inhomogeneous) branching mechanism ip(u,x) of the 
form tp(u,x) — u 2 + k(x)u, where k(x) depends on the density of the obstacles. This work draws 
on similar questions for a single Brownian motion. In the course of the proof, we establish precise 
estimates for integrals of functions over the Wiener sausage, which are of independent interest. 
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1 Introduction 

1.1 Superprocesses in random media 

The purpose of this article is to investigate the behaviour of super-Brownian motion among random 
' obstacles, when the density of these obstacles grows to infinity but their diameter shrinks to zero in an 

appropriate manner. More precisely, let us fix d > 2 and a domain D of K d , and let c : M. d — > [0, oo) be a 
bounded measurable function. For every e € (0, |), let us define an obstacle configuration by 

cn 



T e = |J B{x,e), 



where V e is a Poisson point process on M. d with intensity log(£ -1 )c(x)<iz; if d = 2 and e 2 ~ d c(x)dx if 
d > 3, and B{x, e) denotes the closed ball of radius e centered at x. This Poisson point process is defined 
on a probability space (Jl,^ 7 , P). On a different probability space, let us also consider a superprocess 
^ . {Xt, t £ [0, oo)} with critical branching mechanism ip(u) — u 2 and underlying spatial motion given 

by Brownian motion killed when entering D c U T e . Thus, for each e, the superprocess X s can be seen 
as evolving in a random medium given by T s . A realization of {r e , e S (0,1/2)} will be called an 
environment. 

We wish to understand the behaviour of X e when e tends to zero. As in most works about random 
media, two points of view can be adopted : either we fix an environment (quenched approach), or we 
average over the possible realizations of U e >o T £ (annealed approach). Although the results of this paper 
are set in the quenched framework, the main ingredients of their proofs are "annealed-type" calculations. 
Moreover, the latter approach is also useful in obtaining a better understanding of where the scaling 
comes from and of what the limiting process might be. To simplify the analysis, let us first assume that 
D = R d and let us consider a single Brownian motion £, independent of the obstacles. Denote by P^ the 
probability measure under which £ starts from x. Let us define the random time T e as the entrance time 
of £ into the set r E , that is 

T £ := inf{i > : & G T £ }. 

In addition, for all < s < t, we denote by S £ (s, t) the Wiener sausage of radius e along the time interval 
[s,t], defined as 

S E (s,t)={yeR d : wi \S r -y\<e}= \J (£ r + B(0,e)). 

8<r<t 

r£[s,t] 
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The probability that the Brownian motion £ hits T £ before time t is equal to the probability that the 
centre of one of the obstacles lies in S e (0,t). These centres are given by the Poisson point process V s 
and so, by averaging over the random obstacles and using Fubini's theorem, we obtain 

E[P [Te>i]] =Eo[P[7> £ nS £ (0,i) = 0]] = E \cxp~s d (e) [ c(x)dx], (1) 



where 

log^- 1 ) if d = 2, 



Sd(e) 



z 2 - d if d > 3. 



In the case c = v, the integral in ([I]) is just v times the volume A(5 e (0, f)) of the Wiener sausage, 
whose asymptotics have been well studied owing to their connections with physical problems (see e.g. the 
introduction of [18], |10j or [2]). Note that the large-t asymptotics of \(S E (0, t)) are essentially equivalent 
to its small-e asymptotics thanks to the equality in law : 

X(Sx(0,t)) ®t d / 2 \(S t - y2 (0,l)). 
A classical result of Kesten, Spitzer and Whitman (cf. [8], p. 253) states that, if d > 3, 

lima d (e)A(5 e (0,t)) = k d t a.s, (2) 

where k d = (d — 2)7r d / 2 /T(d/2) (k^ = 2tt) is the Newtonian capacity of the unit ball. The Kesten- 
Spitzer- Whitman convergence result was in fact stated for the large-time asymptotics of \(S e (0,t)), but 
a scaling argument gives the previous statement, at least in the sense of convergence in probability. The 
convergence in ^ also holds if d — 2 (see [H]), with = n. 

It is not hard to deduce from the preceding result that, at least when the function c is continuous, 

limsd(e) / c(y)dy = k d / c(£ s )ds a.s. (3) 



Se(0,t) 



It then follows from (!]) that 

limE[P [T £ >*]] =E [exp-fc d / c(£ s )ds 
L Jo 

This argument, which is due to Kac [9], can be interpreted in the following way. When £ tends to zero, 
the obstacles become dense in R d (at least if the function c is everywhere positive), and the Brownian 
motion £t gets absorbed in the obstacles at rate k d c(t;t)- 

Back to our initial problem about killed superprocesses, the result for a single Brownian particle 
suggests that the sequence X e should converge to the superprocess X* with branching mechanism 
i/j(u,x) = u 2 + k d c(x)u and underlying spatial motion given by Brownian motion. We shall estab- 
lish in this work that the distribution of X* is, indeed, the limit of the distribution of X £ as e tends 
to 0, in P-probability. Here, the distribution of X e is a probability measure on the Skorokhod space 
D Mr (R< i )([0) °°)) °f an cadlag paths with values in Mf(R d ) (the space of all finite measures on R d ) and 
the preceding limit is in the sense of weak convergence. A stronger statement can be made, but only for 
subsequences: if the sequence e n decreases to fast enough, 

X — > X as n — > co, P— a.s. 

Here, ^5 denotes convergence in distribution. Let us emphasize the meaning of this result : except for a 
set of zero P-measure, if we fix an environment, then the sequence of superprocesses X £n evolving among 
these fixed obstacles converges in law to X*. Theorem Q] and Corollary Q] are stated in a more general 
setting, allowing the superprocesses to live only within a domain D of M. d . 

The question we address in this paper was motivated by analogous works on Brownian motion. An 
extensive literature is already available on this topic, reviewed for example in [20]. Owing to the well- 
known properties of Poisson point processes, they seem to be a natural way to encode traps and have 
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been frequently exploited in investigations of the behaviour of Brownian motion moving among "hard" 
obstacles, where the particle is killed instantaneously when hitting an obstacle as described above, or 
among "soft" obstacles, within which the Brownian particle is killed at a certain rate. Our approach is 
close to ideas developed by Kac in [9], whose probabilistic method differs from the analytic method used 
by Papanicolaou and Varadhan [13] in a similar context. Both derive the convergence in the L 2 (P)-norm 
of the semigroup of Brownian motion among random obstacles when the number of obstacles tends to 
infinity but their diameters tend to (recall that P denotes the probability measure on the space where 
the obstacles are defined). Subsequently, Brownian motion among traps was studied in different settings, 
in particular by Sznitman, who devised the powerful method of enlargement of obstacles (see |20|). 

The problem of super-Brownian motion or branching Brownian motion among random obstacles was 
addressed recently by Englander in [4], [6] and [5], the latter paper dealing with soft obstacles. However, 
Englander considers the supercritical case (instead of critical super-Brownian motion as we do) and keeps 
the sizes of obstacles fixed. Within the obstacles, a particle does not die but branches at a slower rate. 
His interest is in the long-time asymptotics of the process and, in particular, the survival probability and 
the growth rate of the support. His techniques are mostly analytic, in contrast with the probabilistic 
tools of the present work. 

1.2 Statement of the main result 

Let us first introduce some notation and construct the sequence of superprocesses X £ from the historical 
superprocess corresponding to a super-Brownian motion on M. d , independent of the obstacles. We refer 
to [l] for more details on historical superprocesses and their applications. If E is a topological space, 
Mf(E) stands for the space of all finite Borel measures on E. 

The (Brownian) historical superprocess can be defined as follows. Let W be the set of all finite 
continuous paths in M. d , and note that M. d can be viewed as a subset of W by identifying x with the path 
of length zero and initial point x. Then, let £ be the continuous Markov process in W whose transition 
kernel is described as follows: If £o = {w(r), < r < s) € W, the law of £t is the law under P w ( s ) 
of the concatenation of the paths (w(r), < r < s) and (£ r , < r < t). The historical superprocess 
H is defined as the superprocess on W with branching mechanism *p(u) — u 2 and underlying spatial 
motion given by £. Thus, H takes values in M. /(W). The super-Brownian motion X° on R d , starting at 
/i € M. /(K d ), can then be recovered from the historical superprocess starting at n (which is viewed as a 
finite measure on the paths of length zero) through the formula 

(X?J) = [ H t {dw)f(w(t)) 

for all / bounded and measurable and all t > 0. Here, (z/, /) denotes the integral of / against the measure 
v. 

We exploit this correspondence between the historical superprocess and super-Brownian motion fur- 
ther to construct the sequence of killed superprocesses which is of interest in this work. Let E be an open 
subset of M. d , and recall the definition of the obstacle configuration r e . For every e > 0, the superprocess 
{Xf' E , t € [0, oo)} is defined from the historical superprocess H via the formula 

{X e t' E -,f)= / H t (dw)f(w(t))l {Vse[Qit]iW ( s ) eEnr e } , 
Jw 

for all / bounded and measurable, and all t > 0. It is straightforward to verify that X e ' E is itself a super- 
Brownian motion with critical branching mechanism ip(u) = u 2 and underlying spatial motion given by 
Brownian motion killed when entering E c UT £ . Furthermore, Xq is the restriction of /i to E n 

Recall that we defined ki = n and kd = r ^ 2 2 ^ 7r rf / 2 for d > 3. We also introduce another superprocess 
X*' E , with branching mechanism ip(u, x) = u 2 + kdc(x)u and underlying spatial motion given by Brownian 
motion killed when it exits E. 

In practice, E will be either Dora bounded open subset of D. When there is no ambiguity, we shall 
suppress the dependence on E in the notation. We choose a sequence e„ such that J2 n I l°g £ n| -1 < oo if 
d = 2, and J2 n £n \ log e ra | < oo if d > 3. For instance, we may fix a > 1 and set e n = exp(— n") if d = 2 
and e„ = n~ a if d > 3. 
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We will use the following notation. 

• P M is the (quenched) probability measure under which H starts at fi € Mf(M. d ) C Mf(W). By the 
preceding correspondence, each superprocess X £,E then starts under P^ from the restriction of fi 
to Ed r°. It will be convenient to assume that X*' E is also defined under P M and starts from the 
restriction of p. to E. 

• To simplify notation, X^' E will be a shorthand for the killed superprocess with parameter e n , and 
p£T ),B will be its law under P M . Likewise, P^ (resp. P;' B ) will be the law of X e ' E (resp. 
under P^. 

• For all t > and x S R , P*^ will be a probability measure under which a Brownian motion £ on 
R d , independent of the obstacles, starts from x at time t. 

• T B := inf {f > : & G £ c }, T £ := inf{t > : & G T e } and T (n) = T £n . 
We can now state our main result. 

Theorem 1. For every fi € M.f{D), P-a.s. 

¥ (n),D ^ p .,£> af , n ^ 

where the symbol =>■ refers to the weak convergence of probability measures. 

As an immediate corollary, we also have : 

Corollary 1. For every fi € Aif(D), the sequence ¥"^ D converges in P -probability to T , * ,D as e tends to 
zero. In other words, for every S > 0, there exists Eq > swe/i that for all < e < sq, 

P[d(V^ D ,¥^ D ) >S] <S, 

where d is the Prohorov metric on Mi(D_m [0, oo)) (here, M.i(Dj^ f ^£>) [0, oo)) is i/ie space o/ a// 
probability measures on D_\4 f ^[0,oo)). 

The rest of the paper is devoted to the proofs of Theorem [T] and Corollary [Tj In Section 2, we 
prove certain estimates for the rate of convergence in J3]), which are of independent interest. These 
estimates are a key ingredient of the proof of Lemma[2]in Section 3. Then, we fix a bounded open subset 
B of D and prove the almost sure convergence of the distribution of X^' B in two steps. First, we 
show in Section 3 that to each k— tuple {t\, . . . ,tk), there corresponds a set of P-measure zero outside 
which (x£ ) ' B ,...,xg ) ' B ) n >i converges in law to (X*^ , X%' ). Second, we prove in Section 4 that, 
with P-probability 1, the sequence of superprocesses X^' B is tight in -DjU/(£>) [0j oo). In Section 5, we 
complete the proof for a general domain D. Starting with a bounded subset of D is required for technical 
reasons, to ensure the finiteness of certain integrals which appear in the proof. 

2 Some estimates for the Wiener sausage 

Let us define the set B\ as the set of all bounded Borel measurable functions c on R d such that ||c|| < 1, 
where ||c|| denotes the supremum norm of c. We have the following result (we write E x for Eq, x in the 
rest of the section): 

Proposition 1. For every t > 0, there exists a constant C = C(t) such that for every e S (0, \], if d = 2. 

f|loge| / c(y)dy-ir [ c(£ s )ds) 



sup sup E. 

c6Bi x6t 2 



Se(0,t) 



C 

< 

I log e | 



and if d > 3, 



sup sup E 2 

c6Bi xEK d 



(e 2 - d ( c(y)dy-k d f c^dsY] < Ce 2 \ loge| 2 . 
v Js e (o,i) Jo ' 1 
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Remark 1. In the case c = 1, the bounds of Proposition^ follow from the known results for the fluctu- 
ations of the volume of the Wiener sausage J12\/ - However, it does not seem easy to derive Proposition^ 
from the special case c = 1. Note that the latter case suggests that the bound Ce 2 \ loge| 2 could be replaced 
by Ce 2 \ loge| if d = 3 and by Ce 2 if d> 4. These refinements will not be needed in our applications. 

Proof of Proposition [l] for d > 3: To simplify notation, we prove the desired bound only for t = 1. A 
scaling argument then gives the result for any t > 0. Let us set 



h(e) = sup sup E x 



(e 2 - d / c(y)dy-k d / c(Qds) 
v Js e (o,i) Jo ' 



As a first step, let us notice that 

c{y)dy 



S E (0,1) 



Also, 



S e (0,l/2) 



1/2 



c{y)dy 



I c(y)dy- \ 

JSAl/2,1) JSAOA 



(0,l/2)nS e (l/2,l) 



c{y)dy. 



MO, 1/2) 



c(y)dy-k d c(Z s )ds = e 2 ~ d 2- d / 2 



S eV5 (0,l) 



i) 



dz 



where £ s = y2 £ s / 2 for all s > and S e (a,b) is the Wiener sausage associated to £. Since the function 
c(z) = c(^=) also belongs to B\, we obtain that 



(' 



2— d 



S s (0,l/2) 



1/2 \ 2 1 1 

c(&)ds) < -h(eV2). 



Likewise, using the Markov property at time 5 and the preceding argument, we have 

1 



E„ 



Se(l/2,1) 



c(y)dy-k d [ c(Zs)ds) 2 ] < jh(sV2). 

Jl/2 ' 1 4 



On the other hand, we have A(S e (0, l/2)nS e (l/2, 1)) = \(S , £ (0, l/2)nS^(0, 1/2)), where £ = £i/ 2 -t-£i/2 
and = £1/2+4 — £1/2 for every t e [0, 1/2], and S' e (0, 1/2), resp. S"(0, 1/2), denotes the Wiener sausage 
with radius e associated to resp. along the time interval [0, 1/2]. Since £' and £" are independent 
Brownian motions, we can use the following consequence of Corollary 3-2 in [11], and of [12], p. 1012: 
There exists a constant Ki(d) > such that for every eg (0, 1/2] 



E 



i^ie 4 , d = 3 

\(S E (0, 1/2) n S £ (l/2, 1)J J < I K l£ *\ loge| 2 , d = 4 



d> 5. 



Coming back to the definition of h(e), and using the triangle inequality in L 2 , the fact that 



Ex / c(y)dy 

L v JS e (0,l/2)n5 e (l/2,l) 

and the preceding inequalities, we obtain 



< E 



A S E (0,l/2)nS E (l/2,l) 



h(e) < sup sup ^ E x 

cSBi x£R d 



(e 2 ^ / c(y)dy + e 2 ~ d c(y)dy - k d c(£ s )ds) 

K JS e (0,l/2) JS e (l/2,l) JO 7 



1/2 



-E 3 



_4-2rf 



< 



S e (0,l/2)nS e (l/2,l) 
1/2 



c{y)dy 



1/2 



{ (^/i(eV2) + 2u(e)) * * + ^(e)} , 



(4) 
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where ^rf( £ ) = e (resp. e 2 \ logs], resp. e 2 ) if d = 3 (resp. d = 4, resp. d > 5) and 



u(e) = sup sup 













1/2 



c(y)dy-k d 

/2) 

Applying the Markov property at time \, we have 



■(Cs)d S ) (e 



Se(l/2,1) 



z{y)dy-k d / c(£ s )d. 

Jl/2 



— sup sup 



E, 



S e (0,l/2) 



1/2 



c(y)dy-k d c(£ s )ds)u(£,£i/ 2 ) 



where 



_2-d 



S e (0,l/2) 



1/2 

c(y)dy - kd I c(£ s )ds 



We now use the following lemma. 

Lemma 1. There exists a constant K 2 > such that for all z e € (0, ^] and c £ Bi 

\v{e,z)\ <K 2 e. 

We postpone the proof of Lemma [Hand complete the case d > 3 of the Proposition. By Lemma [H 
we have 



\u(s)\ < K 2 e sup sup E 2 



Se(0,l/2) 



1/2 



From (@]), we obtain for every e G (0, ^] 

h{e) < ^h(eV2) + K 2 eh(eV2)^ 2 y /2 + K[Me))' 
Let us set g(e) — e^ 1 h(e) 1 ^ 2 . We thus have for e € (0, |]: 

< ( 5 (eV2) 2 + ^ 2 . 9 (eV2)) 1/2 + J^ e - 1 ^(e). 



c(y)dy-k d ^'"cfe)ds) 2 ] 1/2 < ^ £/i(eV2) 1/2 . 



(3) 



Fix r e (1/4,1/2] and set u n = g(r2~ n / 2 ) for every integer n > 0. Rewriting ((5]) in terms of it„ and 
noting that e~ 1 ipd(£) = 1 if d = 3 and e~ lr ipd(£) = o(l) as e — > if d > 4, we obtain for a constant 
K" > (independent of n) 



in+l < {U 2 n + V2K 2 U„) 1 / 2 + = U n (l 



V2if 2 \ 1 ' 2 



V2^2 

2 



It follows that u n < uq + n(K 2 2 1 / 2 + K'{) for every ?i > 0, from which we can conclude that there exists 
a constant K3 such that for all e G (0, 1/2], 



and thus 



g{e)<K z \\oge\ 
h{e)<K 2 e 2 \\oge\ 2 . 



□ 



Proof of Lemma Q} We may assume that z = 0, and we fix c € B\ (the constant K 2 will not depend 
on c). First, we have 



Er 



1/2 1 

kd I c(£ s )ds = k d I dy c(y) 



^ ds 
o {2ns) d / 2 



exp 



(-fy 



(6) 
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Let us define the random times T £ (y) and L E (y) for all e > and y € M d by 



r e (y) = inf{t>0: |& - y| < e}, 
L e (y) = sup{t>0: \&-y\<e}, 

with the conventions that inf = +00 and sup0 = 0. We thus have 



s e (o,i) 



c{y)dy 



dy c(y)P |r e (y) < - 

L Z 



c(y)P < L e (y) < - + / c(y)P r E (y) < - < L £ (y) 



(7) 



On the one hand, 



We have 



dy c(y)P r £ (y) < - < L £ (y) < dy P r e (y) <-< L £ (y) 



dy P c 



T s (y) < 2 < Ml/) 



— E 
= En 



A^5 £ (0, 1/2) n5 , E (l/2,oo) 
a(s e (0,1/2) nfi£(0,oo) 



where S' e denotes the Wiener sausage associated to a Brownian motion £' independent of £ and also 
started from under Pn. If d = 3, it is easily checked that 



Eq 



(5 £ (0, 1/2) nS' e (0, 00)) = 0(e 2 ) 



(8) 



(use the fact that P [y € S' e {0, 00)] = A A 1, together with the bound (3.d) in [12]). If d > 4, 



En 



A(s £ (0,l/2)n^(0,oo)) = e d E A(5i(0,e- 2 /2) n5i(0,oo)) 



0(£ 4 |loge|), if d = 4, 
0(e d ), if d > 5 



(9) 



by [12], p.1010. 

Let v £tV {dz) denote the equilibrium measure of the ball B(y,e), that is the unique finite measure on 
the sphere dB(y,e) such that for every x with \x — y\ > e, 



P x [r £ (y) < 00] = J v £ . y {dz)G{z - x), 



where G(z) = f °°(27rs) _d / 2 exp(— \z\ 2 /2s)ds = Cd\y\ 2 ~ d is the Green function of d-dimensional Brownian 
motion (cy is a constant depending only on d). By a classical formula of probabilistic potential theory 
(see [15], p. 61-62) we have 



< L £ (y) < - 



1/2 



ds I v £ ^ y {dz) 



[2TT S ) d /' 



■ exp 



(-£)• 



It is well known that v £ . y — kdS d 2 7r Eii/ , where ir £ ^ y denotes the uniform distribution on the sphere of 
radius e centered at y. Recalling ([6]), (7j), (JSj) and ([9]), we can write 



\v{e,z) 



kd / dy c(y) 



l ' 2 ds 



(27Ts) d / 5 



(dz)(e 



-Is/172. 



< k d dy 



!/ 2 ds 



(2ws) d / 2 



-K £ ,y(dz) 



e -\z\ 2 /2s _ e -\y\ 2 /2s 



-0(&(e)) 
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where 0d(e) = e (resp. £ 2 | loge|, resp. e 2 ) if d = 3 (resp. d = 4, resp. d > 5). It follows that 

- 1 / 2 ds 



dy 



|j/|<10e 



< 2 



(2-Ks) d / 2 
dz 



|*|<lle 



^e : y(dz) 

l l 2 ds 



e -\z\ 2 /2s _ e -\v\*l1s 



{2ns) d l 2 



-\z\ 2 /2s < 2 



dz G{z) = 0(e 2 ). 



|*|<lle 



On the other hand, we can find constants C and C such that if \y\ > lOe and \z — y\ = e, 



< c 



M 2 -M 2 



exp 



<c £ M exp (. 



4s 



Thus, with a constant K which may vary from line to line, we have 



dy 



\y\>We 



1/2 



ds 



(2ns) d / 2 



< Ke 

= Ke 

< Ke. 



dy \y\ 



\y\>We 



/^„(^)|exp(-!i)-exp(-!fL) 



dy \y\ 



l-d 



\y\>We 



l/(2|y| 2 ) 



ds' s '~d/2-l e -l/4s' 



Combining the above, the proof of Lemma Q] is complete. 
Proof of Proposition [l] for d = 2: Let us define 



h{e) = sup sup E 

ceBi x£R d 



logej 



c(y)dy -iv c(£ s )ds 



Se(0,l) 



By Corollary 3-2 in [IT], we have 



E 



A(£ e (o,i)nS e (i l))' 



< 



|log£ 



The same technique as in the previous case yields 

ft-(e) < ((-h(eV2) + h(eV2) 1/2 sup sup \v(e, 
^2 ceBi ie» 2 

where 



4 ' 



1/2 



|loge| 



v(e,z) = E 2 



I log e | 



1/2 



(10) 



c{y)dy - TT / c(£ s )ds . 

/ S^O, 1/2) JO J 

We now use the following result, whose rather technical proof is deferred to the Appendix : There exists 
a constant K2 such that, for e S (0, 1/2], 

^2 



sup sup u(e,z) ^ 1 , 



(11) 



1/2 



Hence, if g(e) = | loge|/i(e) 1//2 , we have for e e (0, 1/2], 

\2(log£V2) 2 |logeV2| 

From lfl2|). we can use arguments similar to the case d > 3 to infer that the function g{e) is bounded over 
(0, 1/2]. Thus, there exists a constant K3 such that for all e € (0, 1/2], 



K\ . 



(12) 



h(e)< 



K 3 
|loge| : 



□ 
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3 Almost sure convergence of the finite-dimensional distributions 

of X^ B 



In the following, we fix a bounded open subset B of D and consider only the superprocesses killed outside 
B. We therefore suppress the dependence on B in the notation. In particular, T = T B . 
The following proposition is the first step in the proof of Theorem [TJ 

Proposition 2. (Convergence of the finite- dimensional distributions) Let fi £ Mf(B), p £ N 
and t\ < . . . < t p € [0, oo). Then, under 



(v( n ) 



,x^ ) ) ( ^{x; i ,...,x; p ) 



as n — > oo, on a set of P -probability 1. 



Proof of Proposition [2] : We fix an environment. Let p £ N, < t\ < . . . < t p and ft, ■ ■ ■ ,f p £ 
2?b + (R d ) be measurable, nonnegative and bounded functions. In the following, we shall denote (t\, . . . , t p ) 
by t and (/i, . . . , f p ) by f. 

Let [i £ A4f(B). Following the notation in |3J, we have : 



exp-^(X 4 e 4 ,/i) =cxp-(^,wg), 

2=1 
P 

CXp-^^*,^) =CXp-(pt, Wq), 



where w e = (wf (x); t > 0, x £ B) and w* = t > 0, x £ B) are the unique nonnegative solutions 

to the following integral equations: for all x £ B and t > 0, 



(x) + E fv 



U^+E^ / rfs^(6) 2 I{ S <TAT e} =VE t , x /ifeJI^^TATJ 

LJt J i=i L 

i P r 

K(6) 2 + M6K(6)) I{ S <T } =E Et ^ M&J hu<T} 



i=l 



(13) 



(14) 



where by convention E t a; [/(^ s )] = if s < t. By the standard argument of the proof of the Feynman-Kac 
formula, the integral equation (f!4|) for w* is equivalent to 



w* t (x) + E t . 



ds w* s {t s ) 2 e- k *S°^ du 1 {S<T} ] ^^tJmu)e- kd ^ c(iu)du hu<T} 



(15) 



The equivalence of the two integral equations (JT4J) and (flB")) corresponds to the well-known fact that 
super-Brownian motion with branching mechanism ip(u,x) = u 2 + kdc(x)u can also be constructed as 
the superprocess with branching mechanism rb(u) — u 2 and underlying spatial motion given by Brownian 
motion killed at rate kd,c(x). 

Remark 2. Since w £ and w* are nonnegative, fTS\) and {1$ imply that iy| and w e t are equal to zero 
whenever t > t p (recall that by convention, the right-hand side of $13)) or {1$ is zero when t > t p ). 
Likewise, uf (x) — if x £T e C\ B, for every t > 0. 



By integrating over B the difference between (fl3|) and lfl5|) . we obtain : 



dx \wt(x) - wl(x)\ < dx ^E t)!E fi(£,t i )^{u<T}O i {u<T E } - e 



-k d f t 4 c(£„)dti 



(16) 



dx Et 



rf.7' 



ds I {s <T}(I {s <T e} ~ e- k «f° c ^ du )wl(Z s ) 2 _ 
ds e- k «f t Sc ^ du I {s<T} (wUt s ) 2 - <(6) 2 ) 
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Let us start with the third term in the right-hand side of l|T6|). The functions w e and w* are bounded by 
Cf : = ELi ll/ill. hence bounding I {s<T} by and e"^ // c ««)*« by 1 yields 

ds e- k ^ t 3 <^ du l {s<T} (w;^ s ) 2 - <(U 2 ) 

d s i s (£ s )K(6)-<(&)l 

dz \w*(z) - w e s {z)\p s - t {x, z) 



\ dx 




IB 





< 2C f / dx E tj , 

' B 



2C t 



BxB 



ds 



< 2C f / ds dz \w*(z) - w £ s (z)\. 



(17) 



In the preceding estimates, p r (-, ■) denotes the transition density at time r of d-dimensional Brownian 
motion. The last inequality stems from the observation that j B p s -t(x,z)dx = J B p s -t(z,x)dx < 1. 

We next show that the first two terms of lfl6|l converge towards P-a.s. The key ingredient is the 
following result: 

Lemma 2. Let t\ £ [0, oo) and let f £ £>&+(R d ) be a bounded nonnegative measurable function. Then, 
there exists a constant K — K{c, t±, d) such that, for every t £ [0, oo), x £ B and e £ (0, 1/2), if d = 2 



and if d > 3, 



E 



E 



E, 



E, 



/(eownCWT.}-^ 1 *"^)! 2 ! < #11/11 



I log e\' 



/(6x)I {tl <T}(l {tl <T e} -e- k *ti lc ^ du )Y] < K ll/f £ |log s\. 



The proof of Lemma [2] is postponed until the end of the section. Let us temporarily fix t £ [0, t p ] 
Applying the lemma with e — e n , we obtain for every 5 > and every i £ {1, . . . ,p} 



dx > S 



p 













< 



< 



mu)l { u<T } (hu<T (n)} ~ e- k *fr c ^ du )] I dx) 













E/..C 



A(g) 

/A(£?) 2 A' ||/ 4 || 2 <5- 2 llogenl^if d = 2, 
\A(i3) 2 X ||/ 4 || 2 <5- 2 e n |loge„| if d > 3, 



d.c 



which is summable by our assumptions on (e n )n>i- Hence, by the Borel-Cantelli lemma, 



P-a.s. 





Et,a; 







Mtu)hu<T } (hu<T M }-e- k «^ c ^ du ) 



dx^O 



as n tends to infinity. The first term of 1|16|) is bounded above by a finite sum of such terms, therefore it 
converges to P-a.s, for each fixed t £ [0,t p ]. 
Let us set 



Af,t := (m) £!!x [0,t p ] : / dx V^A^ h{tu)l {u <T}{\u<T {n) } ~ ^ kd ^ c{U 
I Jb i=i 



If \\ denotes the Lebesgue measure on R, we have by Fubini's theorem P (g> \i(Af t ) = J P dt P({w : 
(w,t) G A c f t }) = 0, which gives (i) in the following lemma : 
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Lemma 3. (i) There exists a measurable subset i!f,t of £1, with P(f2f,t) = 0, such that for every 
lu 6 f2\f2f ; t, 



^[Z,M&)h**<n(ku<T M }-e- k '# c{Mdu ) 



dx — > as ra — ► oo 



/or aZZ t > 0, except for t belonging to a Lebesgue null subset Tft.^ o/R+. 

(ii) There exists also a measurable subset Cl[ t of ft, with P(fif i t) = 0, such that for every u> € fi\f2f,t, 



ds / dx 
Jb 



I {S <T} (I {s <t ( „ )} ^ e- k « fo «M*")w<$ t fo)< 



as n — > oo 



/or aZZ i > 0, except on a Lebesgue null subset Tf.t,^ o/R+. ifere, u^™* 1 = w e ™ is ifte function given by 
$13\) corresponding to the superprocess X^ n \ 

(Hi) Finally, for all x € B there exists a negligible set £lf,t,o(x) outside which 



{U<T {n) } 



-k d Jo' 1 c(£„)du 



i=l 



and 



Eo,*[j2Mtu)hu<T}<k 

ds I {s <t}(I {s <t ( „ )} - e-^o'*)*)^^) 2 



converge to as n — > oo. 

Both (m) and (m) can be obtained from Lemma[2]in a way similar to the derivation of («). Note that 
in (m), we may replace the integral over [0, oo) by the integral over [0, t p ] (since wi^ = if r > t p ) and 
that the functions Wr are uniformly bounded by Cf . 

The first term of the right-hand side of (fT6| . with e = e n , converges to as n — > oo provided that 
w *tlfx and i ^ ^t.u, by Lemma [3] (i). For the second term, we have 



dx 



E, 



(18) 



da; 



Eq.x 



ds h 



{ S <T}(I{ S <T ( „)} 



< 



ds / dx 
Jb 



E 



0..r 



I{ S <T } (l {s <T ( , l)} -e- fc ^o^)^) W i" + ) t (6)^ 



which converges to as n — » oo by Lemma [3] (ii), ii lu £ f2f ; t and t £ 7f,t,w 

Finally, for cj € (fVt U &f,t) C and i € (Tf,t,uj U Tf,t,^) C , the first two terms of the right-hand side of 
(flG)) converge to as n — > oo. Recalling (fl7)) . we obtain 



K (n) (x) - w*(x)| < &„(*) + 2C f / ds I dz \w { t ' l, (z) - toJ(«)| 

Jb 

where b n (t) -* as n -> oo provided to and t are as above. Besides, for every t, 

\b n (t)\<2\(B)(l+t p )(C f + Cf). 

Set for every t G [0,t p ], 



,(«)/ 



(19) 



G„(t) := / dx |«£l t (z) - <_ t (a:)| 



Then, 



G n (t) <b n (t p -t)+K f / dsG„( S ) 
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where Kf := 2Cf. By iterating this inequality as in the proof of Gronwall's lemma, we obtain for all 
k > 1, n > 1 and t G [0, t p ] 



fc-2 



i=0 



G n {t) < b n (t p -t) + ^2 K f +1 ds i ds 2 ... ds i+1 b n (t p - s l+1 ) 



JO 



KS / dsi / ds 



fc-2 



ii JO 

t 



ds k G n (s k ) 



o Jo 



< b n (t p -t)+J2 K f +1 I d Sl I ds 2 ... I ds l+1 b n (t p - s i+ i) + \(B)K* +1 ^ 



Fix £ > and let k > 2 be such that \{B)K^ < §. For w <£ (fi f:t U Of, t ), b n (r) converges to 
n — > oo except on a Lebesgue null set of values of r, and thus by dominated convergence 

ds\ I ds 2 ■ ■■ ds l+1 b n (t p - Sj+i) — > 
o Jo Jo 

for every t G [0, £ p ] and i G {0, . . . , fc — 2}. In particular, for such w and for every t G [0, t p ], we have 



as 



dsi / ds 2 ... ds i+1 b n (t p - s i+1 ) < 



2k 



for all 7i sufficiently large. If moreover i is such that t p — t € (^Tf,t, w U 7f,t,wJ ; then we have also 
&n(ip — t) < ^ if n is large. Hence, we have G n (t) < e when n is large. Since e was arbitrary, 
we can conclude that for all u and t as specified above, G n (t) converges to 0. Equivalently: for all 
lu e (fif.t U A f ,t) C and i G (ff, t ,u U Tf )tj(1 ,) c , 



hm / cix |iUj (a:) — = 0. 



(20) 



We next consider the asymptotic behaviour of ^^(x) — Wq(x)\. In the same way as in (fl6|) but now 
without integrating over B, we have for every x G B 



U n \x) - w* (x)\ < \E 0tX [j2fi(^)l {U <T } (l { u<T w} ~ e- k *& c^) du) 



i=l 



(21) 



r r°° 

+ |e , x [ / ds i {s<T} (i {s<T( „ )} - e -^/;*)^>W(y2 



2CfE , x 



d s M£ s )K(6)-^ n) (6)l 



Let us fix a; € -B. By LemmaH (m), there exists a P-negligible set £lf,t,o{x) outside which the first two 
terms in the right-hand side of (|2Tj) converge to 0. Besides, for any S > 0, 



Ei 



0..r 



f 

JO 



dsI B (( a )\w* s (Z s )-wi n \Z s )\\= ds dzp s (x,z)\w* s (z)-w{ n Hz)\ 



< 2C f S 



1 



ds I dz \w*{z) - wi n \z)\, 



(2iT6) d / 2 

using the bound p s (x,z) < (27r<5)~ d / 2 if s > 6. If in addition to G ffif,t U f2f t ^ , then by ((20j) and 
dominated convergence (recall that w* s and vanish for s > t p ), 



ds / dz \w*(z) - wi n \z)\ ->0 
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and so limsup \wq (x) — w'q(x)\ < 2Cf8. Since 5 was arbitrary, it follows that lim \wq (x) — w'q(x)\ = 0. 
To summarize, for all x G B and u G (&t,t U ^f : t U fif : t,o(a;)^ (of P-probability 1), 

lim \w { n) (x) - w* (x)\ = 0. 

n — >oo 

From the latter result, we can obtain the convergence of the finite-dimensional distributions of 
towards the corresponding ones for X* . For all x € B, P[w^ l \x) — > Wglz)] = 1 so by applying once 
again Fubini's theorem, we have 

P— a.s., fi— a.e., (x) — > Wq (x) as n — > oo. (22) 

Since the t// n ' are bounded by Cf , dominated convergence and (|22| give 

exp-(/U, w[, n) (-)) -> exp-(/x, «£(•)). 

Our construction from the historical superprocess makes it obvious that X^ is stochastically bounded 
by X°. It follows that the sequence of the distributions of {(x[™\ . . . , X^), n 6 N} is relatively compact. 
Therefore, if we choose a countable set of p-tuples (/i, . . . , / p ) such that the corresponding family of maps 

(fil, . . . , Hp) i— ► exp — X)i=i(Mii /i) is convergence determining, we obtain that (Xj" , . . . , -Xt ) converges 
in distribution to (X£ t , . . . , X t * ) on a set of P-probability 1 (which a priori depends on (ti, . . . , t p )). This 
completes the proof of Proposition [2l 

□ 



Proof of Lemma [2j The quantity of interest vanishes if t > t\, and so we need only consider the case 
t <t\. In that case, 



E 



(E M [/(COWt} (h tl <T e} - e~ kd X 1 c ^ du )]) 



(23) 



d a 



hti<T>}e~ 



kd it 1 c(Su)du + e -k d Jt 1 (c(U+c(C))du 



where £' is another Brownian motion, independent of £, T" and T' e are defined in an obvious way and 
we have kept the notation Pt,x for the probability measure under which the two Brownian motions start 
from x at time t. Recall that S £ (s, t) denotes the Wiener sausage of radius e along the time interval [s, t] 
associated to £, and define S' e (s,t) in a similar way. Then, 

E[l {tl<Te} I {tl<T , } ] =P[7> £ n(S e (i,i 1 )U^(M 1 )) = 0] =exp\-s d (e) [ c(y)dy\ 

L J Se(t,tl)US' e (t,tt) ' 

and similarly 

E[l{ 4l <T t }] =exp{ - s d {e) / c(y)dy\. 

L Js e (t,tx) J 

Set ti = ti — t. By Fubini's theorem and a simple symmetry argument, the quantity in (|23|) is equal 
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to 



E, 



W}W^/&0/(O{2e- fc ^ Ic( ^ 



<2ii/rE M 



-^( e )/ Se ( t , tl )us^( t , tl ) c fe)^ _ g-^/t' 1 (c(€ u )+c(?:))d«~ 



+ ll/l| 2 E t , :c 


e -^(£)/ Se(t-tl)us , (t , tl) c(y)dy _ e -k d j^ (< 






<2||/|| 2 E ,, 


fcrf / c{£,u)du - s d (e) / c(y)dy 


+ ll/H 2 Eo,, 


sd(e) / c(y)dy 




■/ JS e (0,ti) 




•/S e (0,*l) 



c(y)cfy - s d (e) 



S£(0,tl) 



c(y)dy| 



c(y)dy-k d c(£, u )du-k d c(£' u )du 



ti 



<4\\f\\ 2 E . x | k d I c(£ u )du - s d {e) 

'S e (0,ti) 

+ ||/|| 2 || C || Sd (£)E ,4A(^(0,t 1 )n^(0,i 1 ))], 

where in the second inequality we used the bound \er x — e~ y \ < \x — y\ for x,y > 0. 
On the one hand, by p] (d = 2, 3) and [12J p.1009-1010 (d > 4), we have 

|loge| Eo,x[A(S e (0,ti)n^(0,ti))] =0(\\oge\- 1 ) if d = 2, 

£ - 1 E ,4A(S' £ (0,ti) nS^(0,ti))] =0(e) if d = 3, 

e 2 - d E ,4A(5 e (0,?i) n £(0,ii))] =0(£ 2 |log £ |) if d > 4. 

On the other hand, 



(24) 



(25) 



Eo,x 



k d / c(£ u )du - s d (e) 



/ c(y)dy 


< E ,a; 




/S e (0,tl) 







c(y)dy 



S e (0,ti) 



1/2 



(26) 



Proposition Q] ensures that the right-hand side of (|26|) is bounded by i^||c|| |loge| 1 if d — 2 and by 
K\\c\\ e\ loge| if d > 3. Together with ([24|) and |25|) . this completes the proof of LemmaO 



□ 



4 Tightness of the sequence X^ ,B 

Let C 2 _(IR' 1 ) denote the set of all nonnegative functions of class C 2 on K d . By Theorem II. 4.1 in [14], 
the tightness of the sequence of the laws of the superprocesses X^' 3 will follow if we can prove that 
the sequence of the laws of (X^' 3 , 4>) is tight, for every cj) g C^(R d ) with compact support. Note that 
condition (i) in Theorem II. 4.1 of [H] holds thanks to the domination X^' 3 < X°. Recall that A" is 
the usual super-Brownian motion without obstacles. 

Let us first introduce the P-negligible set 9 C outside which the desired tightness will hold. 

Definition 1. (Good environments) Let O be the union over all choices of the integer p > 1 and of 
the rational numbers qx, . . . , q p of the P-negligible sets on which the sequence (Xg\...,x£>) does not 
converge in distribution to (X* t , . . . , X* p ) as n — > oo . We call good environment any environment which 
does not belong to 8. 

To simplify notation, we again write X^ for X^ n ^ B (as in the last section, B is fixed) and prove 
tightness only on the time interval [0, 1]. Let us fix g C\ (R d ) with compact support. The tightness of 
the sequence {X^ n \ 4>) is a consequence of the following lemma. 
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Lemma 4. If uj £ 6, then for every e > there exists k — k{e) > 1 and uq = no(uj,e, k) such that for 
all n > n , 



k-l 



[U{ sup \(X<; n \<j>)-(xf,<j>)\>e} 
' ■• =o i<*<— 



< e. 



Lemma [4] easily implies that the sequence (X^ n ',(/)) is tight. Indeed, let us fix a good environment 
and Tj > 0. By LemmalU there exist k{rj) and no(ui,i], k) such that for all n > no, 



fc-i 



,[U{ sup i \{Xl n \ct>)-{xf^)\>l} 



< V- 



(27) 



If n > n is fixed, on the complement of the event considered in lj27j) . we have for every s,t £ [0, 1] 



t-s < 



(xl n \<f>)-{XW,cj>)\<r, 



(n) 



and therefore w({X^, 0), r, 1) < n, using the notation of Ethier and Kurtz [7] for the modulus of 
continuity of the process (X^ n \cf>), Thus, for all n > no, 



k— 1 

i 1) < „] > P„[( (J { 4 sup |<X^,0) - <X<j»>, 0>| > |}) 



o 



> 1 - 77. 



In addition, is bounded so the first condition of Theorem 3.7.2 in [7] is trivially fulfilled, hence Corollary 
3.7.4 of [7] implies that for any good environment, the sequence of the laws of (X^ n \ 0) under P M is tight. 
Let us now turn to the proof of LemmalU 

Proof of Lemma [4l We fix a good environment. Let e > 0. The process ((Xf,(j)))t>o is continuous, 
therefore there exists ko(e) such that for all k > ko, 



sup 

0<i<fe-l 



(XUA)-{X\ 



> 



< 



There exists K = K(e) > 1 such that 



sup (X°,1)>K 

0<t<l 



< 



(28) 



(29) 



By a trivial domination argument, the bound l(29|) remains valid if we replace X° by X^ (in fact for 
any environment). In the following, we fix the constant K > 1 such that ([29)1 holds. 
We now have the following result : 

Lemma 5. There exists a constant C = C{<p 1 K) such that for every integer k > 1 and every measure 
7 G M f (R d ) satisfying (7, 1) < K, 



sup 

0< S <-1 



(x° s ,<i>)-{x°,, 



£' 

> 2 



< 



k 2 



The proof of Lemma [5] is deferred to the end of the section. Let us define 



A n = { sup (Xt\l) >K] 
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Then, 



k-l 

MU{. 



sup |{XW,^)-<XW, 

<t<*+i " 



> e 



sup 

0<i</c-l 



(*5U>-<4 n) 

+ pjA^n{ sup { sup ((lf>MX<^)))> E ) 

— — fc — — fc 

+ F,kn{ sup { sup ({iW,^-(xW^>)}>4n{ 

fc — '•— fc 

From (|29)l . we have 

On<|- 

Moreover, from the definition of a good environment, if fc > fco, 



> 



sup 

0<i<fc-l 



{X%\,<j>)-(X 



(n) 



< 



limsup& n < 



sup 

0<i<fe-l 



, 4>) - (xj , ( 



£' 

> - 
- 2 



£ 

<3' 



by (f28|) - Thus if fc > fc , there exists n (w,£, fc) such that for all n > n , b n (k) < |, Then, 



A— 1 



<]Tp4<X ( / i) ,1)<^; sup (<X< n \0>-<Af\0>)> e 



i=0 
k-l 



: <t< 



EM 1 



i=0 



{<x<"\i><jr} F xf 



sup ((x ( ^,0)-<xr»)> 



(n) 



0<t<, 



The last equality is obtained by applying the Markov property to X^ at time i. By a domination 
argument, we have for all 7 S A4/(R d ) such that (7, 1) <K, 



sup - (ir,*) > e 

0<t<i 



< 



sup «X t V)-(X °,<£)) >e 

0<t<i 



< 



fc 2 



by Lemma [H It follows that 



, C C s 



if fc > fci(e). Finally, 

fc-i 

rf n < V [ sup (X t (n) , 1) < K- sup ((X 

(X%\,<I>)-{X 

k k 

We fix i E {0, . . . , fe — 1} and consider the stopping time 

T,:=inf{t>^: (X t (n) » < (xf, 0) - e} 



;<t<- n 

sup 

0<z<fc-l 



;<*< — 
-(«) 



£ 

< - 

- 2 



Then, the z-th term of the previous sum is bounded by 



i + 1 



< ^ , (X<I l) , 1) < JC, sup ( (X f w , 0) - (X 



,(„) 



> 



E„ 



{T 4 <^1, <X^ n) ,l><K} -4" 



sup ((xw,^-(xr,0))> 



(n) 



0<t<i 
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by the strong Markov property at time 2$. Using Lemma [5] once again, we see that this quantity is 
bounded by C/fc 2 , hence for k > fci(e), 

C e 

Combining the preceding estimates, we obtain that for k = fc (e) V fci(e), and every n > no(w,e, fc), 



?J|J{ ™p \{x[ n \<t>)-{x 



(n) 



< e. 



This completes the proof of Lemma [H 



□ 



Proof of Lemma U Let 7 G .M/(M d ) be such that I7I := (7,1) < K. Recall that the process 
((X t °,l)) t >o is a martingale. From the maximal inequality applied to the nonnegative submartingale 



sup (X t °, 1) > 2K 
o<t<l 



< 



1 



We claim that 



E„ 



((X|,l)-| 7 |) 



(2^~|7l) 4 7 

24, 12 



«x?,i)-i 7 ir 



k 3 



(30) 



To prove this claim, recall that Y t = (X®, 1) is a Feller diffusion, whose semigroup Laplace transform is 
given by 

E[exp —XY t \ Y =y]= exp ( - -^L^J , 



for A > 0. Thus, 



E[exp-A(y t - y)| Y = y] = exp (^) = 1 + Xhy X 3 t 2 y + Xh^y + + o(A 4 ), 

as A —>■ 0. From this expansion of the Laplace transform, we derive that 



E 



(Y t -yf\Y = y 



2At i y + 12t 2 y 2 , 



which proves our claim f30|) . It follows that 



sup (X°, 1) > 2K 
o<t<i 



< 



i2MCM±2) 

(2K- | 7 |) 4 fc 2 



Let us denote by Ax,k the event { sup 0<t< i (X t °, 1) > 2K} and by Bx,k the event { sup 0<t< i | (X®, 
4>)\ >§}. Then, 



,[B K ,k] < ^[A K ,k] + ^i A K,k n Bjf, fc ] < p + P 7 [^,fe n Sjr, fc ] 



(31) 



where cq is a constant depending on K. 

In addition, M t := (X^,<j>) — (Xq,(/)) — J Q dr (X°, ^A0) is a continuous martingale with quadratic 
variation 2 J Q dr (X°,0 2 ). By the Dubins-Schwarz theorem (see Theorem V.1.7 in [E]), there exists a 
standard one-dimensional Brownian motion W such that M t = W / (A/) t for all £ > a.s. On the event 



A C K k , we have 



and 



* W, -A#)| < f||A<j||il- < ^ if f e [O.fr 1 ] 



17 



where c\ and C2 are constants depending on <f> and on K. Choose k such that c%k 1 < \e for every 
k > ko. Then for all k > ko, 

r f* 1 e 

P 7 [A^ fc n fl^fe] = P 7 n { sup |M t + / dr (X r °, -A0) I >-} 



0<t<l/fe 



<] 



7 



A^ jfe n{ sup |M t |>-} 



0<t<l/fc 



<P 

C3 



sup |W f | > - 

0<t<(c 2 /fc) 4 - 



< 



k 2 ' 



where c 3 is a constant depending on (j>, K and e. Together with lj3lj) . this completes the proof of Lemma 

m 



□ 



5 Proof of Theorem GQ and Corollary [T] 

The proof of Theorem [TJ in the case when D is bounded is easy from the results of the previous sections. 
Let us take B = D and let £ denote the union of the P-negligible set on which there exist rational 
numbers 1 1 , . . . , t p such that (X^' , . . . , ) does not converge to , . . . , X^ ) and of the P-negligible 
set on which the sequence X^ n ^ B is not tight. The set £ is also P-negligible and on £ c , Theorem 3.7.8 

of [7] allows us to conclude that X^ n ' — * X* when n — * oo. 

We can now use the previous result to complete the proof of Theorem Q] when D is a domain of M. d 
which is not necessarily bounded. 

Proof of Theorem [T] for a general domain D: Let /i be a finite measure on D and suppose first 
that the support of (i is bounded. Under P^,, the superprocesses l'"' are stochastically dominated by 
the superprocess X°, whose range 

K(X°) = (J suppX t ° 

t> o 

is almost surely compact since its initial value has compact support. Consequently, for every e > 0, there 
exists a bounded open subset B of D containing the support of /z such that, for every environment and 
every n > 1, 

and 

V^[Il{X*' D ) c B] > l-e. 
From these inequalities, we can deduce that 

d(F^' D ,F^< B )<2e, n>l; (32) 
d(P^,P^ B )<2 £ , (33) 

where d is the Prohorov metric on M.\{D Mf ^ D ) [0, oo)). By the results of the last two sections, with 
P-probability 1 there exists an integer no(uS) such that for all n > no, 

d(F^' B ,Fl' B ) < e. 

Together with p^ ]) and ([35 ]l . this yields 

^(pW.- ^*. ) < 5e for all n > n (uj), 

hence we can conclude that F^' D converges towards P* :D on a set of P-probability 1. 
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Finally, if the support of ix is unbounded, we can replace lx by the measure fx defined as the restriction 
of (X to a large ball centered at the origin. Using once again the domination of X^' D (for all n > 1) and 
of X*' D by X°, the law of X^' D under P M can be approximated uniformly in n by the law of X^ n '' D 
under P/j, and similarly for X* ,D . The desired result then follows from the bounded support case. We 
leave the details to the reader. 

□ 

We end this section with the proof of Corollary [TJ 

Proof of Corollary [TJ Let us argue by contradiction and suppose that there exist 5 > and a sequence 
{e/c, fceN} decreasing to zero such that for all k > 1, 

P[d(F?' D ,F;> D ) >S] >5. (34) 

By extracting a subsequence, we can always choose eu such that 

oo 

^llogefel- 1 <oo if d=2, 
fc=i 

or 

oo 

y^gfc | log£fc| < oo if d > 3. 
fe=l 

But the latter condition is the only requirement for the sequence of superprocesses X £k ' D to converge in 
distribution to X*' with P-probability 1, yielding a contradiction with (|34|) . 

□ 

A Appendix: Proof of (1111) 

The bound ifTTj) is a consequence of the following lemma. 

Lemma 6. There exists a function ip : M. 2 — > [0, oo] such that J R2 dy ip(y) < oo and for every y £ M 2 and 

ee(0,i), 

P [y 6 5,(0,1)]--^- / dspAy) <jr%, 
|loge|7o |loge| 2 

where p s (y) = (27rs) _1 exp{ — |?/| 2 /(2s)} is the Brownian transition density. 

Remark 3. The convergence of |loge|Po[y G 5 e (0, 1)] towards ir J Q ds p s (y) as e tends to was first 
obtained by Spitzer \17$ . See also J 11)/ and \19$ for related results. 

Before proving Lemma [SJ let us use it to derive ifTTj) . If c is a bounded nonnegative measurable 
function on IR 2 such that ||c|| < 1, then for every e € (0, =] and z € R 2 , 



E, 



loge| / c(y)dy-TT [ c(^ s )ds] = [ dy c(z + y) (\ log£|P [y G S £ {0, 1)] - tt / ds p s (y) 

JS e (P,l) JO J JO 

< |loge| _1 / dy <p(y), 



which is the desired result. 

Let us hence establish Lemma[H The following proof is due to J.-F. Le Gall (personal communication). 

Proof of Lemma [6} If \y\ < 10e, simple estimates show that the bound of Lemma [6] holds with 
ip(y) — C((log|y|) + l) for a suitable constant C. So we assume that \y\ > lOe. We put 



(y)=mi{t>0:\Zt-y\<e} 
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in such a way that {y £ S e (0, 1)} = {r e (y) < 1}. Let a e be an arbitrary point of the circle of radius e 
centered at the origin, and 



ds I 



{\U<e} 



A straightforward calculation gives 



f(e)=e 2 \\oge\+0(e 2 



(35) 



as e — > 0. 

Lower bound. An application of the strong Markov property at time r e (y) shows for every u E (0, 1], that 



Er 



ds 



■{\is-v\<e} 



On the other hand, 



and thus 



E, 



ds 



-{\£s-y\<e} 



<Po[re(y)<u]f(e). 

f 

ds / dz p s (z), 

J\z-y\<e 



E, 



ds 



< 



{|6-tfl<s} 

ds 
2ns 



- ire 2 / ds p s (y) 



\z-y\<e 



ds 



o 

e: 3 



\*-v\< 
Jo ^ eXP t-^/ 



\z\ 2 -\y\ 2 



2s 



<e 3 



(36) 



where the function 



= 1 2/1 



A 1 ^ r |y| 2 i 

X ^ exp t-17/ 



(37) 



is easily seen to be integrable over R 2 . 

By combining the preceding estimates, we arrive at 

7T£ 2 f u e 3 

P [rM <u]>—J^ ds Ps ( y ) - y^y*i(y) 

and using lj35|) it readily follows that 

P o [yeS e (0,l)] -jt^-t / ds Ps (y)>P^ 
|loge|7 |loge| 2 

with a nonnegative function ip\ such that J dy ipi(y) < oo. 

Upper bound. This part is a little more delicate. We rely on the same idea of applying the strong Markov 
property at time r £ (y), but we need to be more careful in our estimates. For every v > 0, we have 

= E n I{ Tefe) <i +1 , } E 5TE(H) / dr I 



En 



ds I 



{\ts-v\<e} 



- En 



l {r c (y)<l+v} 





rl+v-T E (y) 

/ dr P as [\£ r \ < e] 

Jo 



{\Sr-y\<e} 



s=l+v-T e (y) 



where a £ is as previously a fixed point of the circle of radius e centered at the origin. We can rewrite the 
previous expression as 



Er 



i-X+v "I rl+v 

/ dr I {Te (y)<i+v-r}Va s [|fr| < e] =/ dr P Q [r e (y) < 1 + v - r] P Qe [|f P | < e] 
Jo J Jo 



20 



We apply this calculation with v = v e = | loge| 1 . For r £ [0, v s ], Po[T e (y) < 1 + v £ — r] is bounded 
from below by Pq [T e (y) < l] , and thus 



Po [t e (j/)<1] / dr P ae [|£ r | < e] 
Jo 



<E: 



[ {\Hs-y\<c} 



dr PaMA <£]Po[r e (2/) < 1 + v e - r] . 



From the bound (|36|) . we have 



Eo 



() 



l {\is-y\<e} 



<ne z dsp s (y) +e 6 ^ 1 (y) +v E e 2 e 
'o 



2 p -|y| 2 /io 



On the other hand, by (|37|) . 

"l + l> e 

dr Pa E [|£ r | < e] Po [t e (y) < 1 + v e - r] 



>/ 'dr P ae [|£ r | < e ]-£_(,r / ' dsp.(tf) -e*i(l/) 



7TS 



/(e) > Jo 

(j[ " E ^P ae [|C,| <e])(^ dap a (i/)-~¥i(y) 
dr P«.[|fr| <e] / dsp s (y) 



2 . 



— r 



Straightforward estimates give 

-l+v e 



and 



where 



J ° drP ae [\£ r \<s] - £ 2 (ilog|log £ | + 0(l) 



1 + V E ,1 

drP ae [\&\ <e] / dsp s (zj) <e 2 * 2 (2/), 

e t/l+u e — r 



*a(») = ^ ds log {jzt;)p>(v) 
is integrable over R 2 . Summarizing, we have 



Po[re(y) < 1] J dr P Qe [|£ r | < e] < (ire 2 J ds Ps (yfj . y . 



i (|logllog £ ]-Jf) £ ; 



I log £■ | 



Finally, it is easy to verify that 



//(:')= / drP,,.[|;,| ; >£ 2 (|log£|-ilog|loge|-^'). 
and so by dividing the two sides of the previous inequality by g(e), we obtain 



Po[r e (»)<l] f 



ds p s {y) + 



¥>2(y) 



|loge|Jo |loge| 2 ' 
with a function y>2 such that J <p2(y)dy < oo. This completes the proof of Lemma [6J 



□ 
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